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SECTION  1 
INTRODUCTION 


Interactions  between  the  plasmas  in  space  ami  the  surface  ami  various  subsystems 

of  spacecraft  are  very  complicated  and  have  been  the  subject  of  considerable  study  over 

the  past  several  years,  f-lectrostatic  charging/1 ,2’^  for  example,  of  a  spacecraft's 

surface  can  result  in  discharges  which  can  cause  electromagnetic  interference, 

degradation  of  surface  materials,  and  failures  of  sensitive  components.  Techniques  to 

influence  plasma-spacecraft  interactions,  such  as  on-board  plasma  generators  and 

(3) 

conductive  coatings  for  dielectrics,  are  also  being  actively  studied. 

The  plasma  environment  of  space  can  be  partially  simulated  in  the  laboratory 
using  low-temperature  plasma  generators  for  studies  of  phenomena  in  the  ionosphere  and 
low-earth-orbit  or  combinations  of  electron  and  ion  beams  to  simulate  the  conditions  in 
high-altitude  orbits.  Several  small  space  plasma  simulation  laboratories^'^’’  111  and  a 
few  large-scale  facilities  are  in  operation^51  or  being  planned.^  Laboratory  simulation, 
however,  is  necessarily  only  a  partial  re-creation  of  the  actual  environment  to  which  a 
spacecraft  is  subjected. 

The  selection  of  the  plasma  generators  or  beams  to  simulate  the  space 
environment  is  now  based  on  intuitive  as  well  as  scientific,  engineering,  and  economic 
grounds.  The  simulation  often  represents  only  the  most  extreme  case  expected  for  a 
given  spacecraft  component.  There  are  presently  no  established  techniques  for  selecting 
a  laboratory  plasma  environment  to  simulate  the  measured  or  postulated  properties  of 
plasmas  in  space. 

The  object  of  this  work  is  to  investigate  some  mathematical  techniques  which 
could  be  used  to  choose  the  parameters  of  monocncrgotio  beams  to  simulate  space 
plasmas.  The  moderate  temperature  plasmas  of  geomagnetic  substorms  serve  as 
examples  for  simulation,  since  they  arc  known  to  cause  electrostatic  charging  on 
geosynchronous  satellites.  The  multikiloolcetronvolt  energies  and  densities  of  a  few 
particles  per  cubic  centimeter  require  their  simulation  by  monocncrgetic  beams  rather 
than  by  low-energy  plasmas  with  a  continuous  energy  spectrum. 
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The  plasma  environment  of  space  Is  characterized  by  a  wide  variety  of  particle 
energies,  fluxes,  species,  and  spectral  shapes.  The  particle  spectra  vary  with  position  in 
space,  time,  and  solar  activity.  Models  of  the  environment  have  been  developed  in 
various  degrees  of  complexity,  ranging  from  the  definition  of  average  plasma  properties 
such  as  density  and  temperature  at  a  given  altitude  to  presentations  of  detailed  spectra 
of  "typical"  plasma  injection  events  recorded  by  instrumented  satellites. 

A  space  simulation  facility  is  necessarily  limited  by  engineering  considerations  to 
providing  a  few  charged  particle  beams  to  simulate  spacecraft-environment 
interactions.  Presently,  the  parameters  of  the  electron  and  ion  sources  are  selected  to 
provide  only  a  rough  simulation  of  the  plasma  environment.  The  energy  and  current 
densities  provided  arc  often  chosen  to  represent  the  most  extreme  case  envisioned  for  a 
given  spacecraft  or  component.  Apparently,  no  quantitative  techniques  now  exist  to 
measure  the  "quality"  of  a  beam  simulation;  and  no  methodology  has  been  developed  to 
specify  the  number,  relative  current  densities,  and  relative  energies  of  a  set  of  charged 
particle  beams  designed  to  simulate  a  given  plasma  spectrum. 

In  this  section  we  examine  techniques  which  can  be  used  to  specify  the  parameters 
of  multiple'  monoenergetic  charged  particle  beams  which  would  provide  a  mathematically 
correct  and  physically  plausible  simulation  of  a  given  plasma  environment.  The 
techniques  are  based  on  the  piecewise  reproduction  of  the  shape  of  distributed  energy 
spectra  or  by  matching  various  averages  of  the  velocity  distribution  functions  by  the 
monoenergetic  beams. 

In  this  study  we  assume  that  the  space  plasma  to  be  simulated  is  of  high  enough 
energy  and  low  enough  density  so  that  collective  effects  in  the  plasma  can  be  neglected. 
More  precisely,  the  !>ebyc  length  of  the  plasma  is  considerably  greater  than  typical 
dimensions  of  a  spacecraft.  This  assumption  is  justified  for  the  space  environment 
outside  the  plasmasphere  during  geomagnetic  substorms  when  strong  spacecraft  charging 
events  are  recorded. 
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2.2  PIECEWISE  SPECTRAL  REPRODUCTION 


1 


The  simplest  and  most  obvious  method  to  simulate  a  distributed  spectrum  is  to 
break  the  spectrum  into  several  bands  and  provide  monoenorgetic  beams  with  appropriate 
currents  and  energies  to  reproduce  the  distribution  in  a  "piecewise"  manner.  A  very 
close  reproduction  of  the  distributed  spectrum  can  be  made  in  this  way,  provided  there  is 
a  sufficient  number  of  available  beams. 

With  a  limited  number  of  beams,  a  problem  arises  on  the  choice  of  the  energy 
boundaries  between  the  parts  of  the  spectrum  to  be  simulated.  Possible  choices  include 
fractions  or  multiples  of  the  average  energy  (temperature)  or  velocity,  or  boundaries 
which  divide  the  particle  flux  into  equal  fractions  of  the  total  flux.  A  given  spectrum 
may  also  be  divided  to  account  for  particular  feutures,  such  as  a  high  energy  "tail"  of  the 
distribution  function. 

The  principles  involved  in  piecewise  spectral  reproduction  can  be  illustrated  by 
considering  a  Maxwellian  distribution  of  particle  energies, 


f(R) 


2n 

/7 


'(kT 


-  3/ It 


(2-1) 


where  n  is  the  number  density,  kT  is  the  temperature,  and  E  is  the  kinetic  energy  of  the 
particles. 

The  differential  energy  spectrum  of  current  density  crossing  an  arbitrary  surface 
is  given  by 


<it 

dt: 


exp 
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where 
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4  \  ii  m  / 
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is  the  total  current  density;  q  and  m  are  the  charge  and  mass  of  the  particles. 

Integrating  Eq.  (2-2)  over  a  range  of  energy  bounded  by  F'. (  and  H,^,  we  find 
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This  current  density  must  he  supplied  by  u  monoenergetic  beam  with  an  energy 
between  E  and  H9  to  simulate  the  corresponding  p.art  of  the  distributed  spectrum.  The 
energy  of  the  beam  can  be  chosen  in  a  number  of  ways;  a  relatively  simple  choice  is  to 
use  the  value  found  by  averaging  over  the  differential  energy  spectrum  of  the  current 


density. 
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Inherit  ion  of  Kq.  (2-5)  gives. 
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Table  2-1  gives  values  for  i(Ej,li.,)/j  and  EOv^Kg)  for  the  ease  of  a  10  keV 
Maxwellian  spectrum  divided  into  four  ranges  of  energy  with  boundaries  at  0,  7.5,  15,  and 
.10  keV. 


tabff  2-i.  imi; riiwisi:  reproduction  of  maxwfffian 

SPECTRUM  BY  FOUR  BEAMS 


Maxwellian  Temperature  -  10  keV 


Fni-rgy  Boundaries 
i : I ,  Iv>  (keV) 


Normali/.od 
Current  Density 

j(Ki.i-:2)/io 


Beam  Energy 
l'.(F,  ,F2)  (keV) 


0,  7 .5 

0.171 

4.0  8: 

7.5,  15 

0.2(i‘« 

11.20 

15.  10 

0.359 

21.19 

11) 

0.199 

42.53 
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2.3  MOMENT-MATCHING  TECHNIQUES 
2.3.1  Velocity  Moments 

A  plasma  can  oe  characterized  by  various  averages  of  the  velocity  distributions  of 
its  constituent  particles.  In  general,  the  "velocity  moments"  of  a  given  <:i-  ••ibution 
function,  f(v).  are  defined  by 


<r> 

/  .<• 


f  (  V  )  V  '  ll  V 


k  -  0, 1,2, - . 

9 

where  the  4  7rv“dv  term  represents  an  infinitesimal  element  in  (isotropic)  velocity  space. 

The  velocity  moments,  M^.,  can  be  related  to  physical  averages  for  several  values 
of  k.  For  example,  M^,  M^,  M9,  and  M ^  are  related,  respectively,  to  the  average  number 
d  nsity  sN>,  particle  flux,  <NF>,  pressure,  <-P  and  energy  flux,  vEE^,  of  the  given 
particle  type  in  the  plasma. 

>'  =  <  N  >  =  n  (3-3) 


/S  vr\  1/: 

Mj  =  4t  <NF>  -  n  <v>  =  n  (”  ) 
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2  in 
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r  >  =  -,r  11  ■  ) 

b  V  ii  ta  / 


(2-111) 


i  =  8:1 

‘3  m 


„  /;:n,  1,: 

<0"  “  2  "(,  ,,  ) 


(2-11) 


M  _  1 5  2  /  8  kT  \ 

M4  -  64  n  11  1 7  m  / 


(2-12) 


5/  ’ 

,  3  2  /  8  kT  \  ' 

f  '  r  ’i  a  ( -  ) 

J  o  \  n  m  / 


(2-13) 


-11- 


I'lu'  ;ivh'iij;('  speed,  •  v  •,  in  liquation  2-‘J  is  defined  by 


<v>  '  Mq  12-11) 

The  expressions  on  the  right-hand  side  of  l.quations  (2-8)  -  (2-13)  are  given  for  the  case 
of  a  Maxwellian  velocity  distribution, 


i  (v) 


(2-15) 


where  n.  :n,  and  1'  are  respectively  the  number  density,  mass,  and  temperature  of  the 
particles  and  k  is  Moit/.mann's  constant. 

Average  ami  1 1  MS  "Temperatures" 

A  useful  method  for  characterizing  a  non-Maxwellian  plasma  is  to  define  effective 
temperatures  whieli  are  related  to  ratios  of  the  velocity  moments/'^  The  average  and 
KMS  temperatures  are  given  by 
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The  two  temperatures  are  equal  when  the  velocity  distribution  is  Maxwellian. 

2.  ?  Monoonergel  ie  beams  to  Match  Velocity  Moments 

\  techni(|iic  to  simulate  a  plasma  with  a  distributed  velocity  distribution  is  to 
ehoosc  the  velocities  and  particle  densities  of  mononergetie  beams  so  that  their  velocity 
moments  match  those  of  the  plasma.  Under  these  conditions,  the  average  parameters  of 
the  beams,  such  as  number  density,  pressure,  or  energy  flux,  are  equal  to  those  of  the 
plasma  component  under  simulation. 
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In  general,  a  single  beam  can  match  two  moments  of  the  distributed  spectrum,  so 
that  two  beams  can  match  four  moments,  three  beams,  six  moments,  etc.  As  discussed 
in  Section  2.3.4,  it  is  also  possible  to  overspecify  the  problem  and  use  more  than  the 
minimum  number  of  seams  to  match  a  given  number  of  velocity  moments. 


Single  Hearn  i'nergv 


A  monoenergetie  beam  ean  match  two  moments  according  to  the  simultaneous 
equations, 


n 


!>'!> 


M. 
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where  n.  and  v.  are  the  density  and  velocity  of  the  beam  particles. 

For  example,  when  the  zeroth  (number  density)  and  second  (pressure)  moments  are 

chosen, 


n.  =  n 
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(2-19) 


or,  in  terms,  of  beam  energy,  F.^, 


F'b  2  k  TAV 


(2-20) 


If  the  first  (number  flux)  and  third  (energy  flux)  moments  are  used, 
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For  ft  Maxwellian  plasma,  where  T^y  TRMS  =■  T,  l-cj.  (2-24)  simplifies  somewhat, 
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The  beam  densities  and  energies  are  then  found  to  be 


n  =0. 2.82  n 
=  0.618  n 


(2-27) 


E  =  3.007  T 
F-2  =  0.568  T 


(2-28) 


Three  Beam  Energies 

Six  moments  of  the  distributed  spectrum  can  be  used  to  compute  the  densities  and 
velocities  of  three  monoenergetic  beams.  No  analytical  solutions  have  been  found  for 
this  case,  but  iterative  techniques  can  be  used  to  find  solutions  of  the  set  of  six 
simultaneous,  nonlinear  equations. 

As  discussed  in  Appendix  B,  the  beam  velocities  and  densities  can  be  found  in 
terms  of  the  average  speed  and  density  of  the  plasma  particles.  For  the  case  of  a 
Maxwellian  plasma  with  temperature,  T,  the  beam  densities  and  energies  are 


=  [0.087,  0.588,  0. 325 J  n 


1,2,3 


El,2,3  =  t4-93**  1-65/,  0.303]  T 


(2-21) ) 
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I  >i ffiTi-nt  values  will  be  found  for  other  types  of  velocity  distribution  functions, 
Intt  t ho  method  used  to  compute  the  Maxwellian  results  is  general  for  ull  realisti" 
spectral  shapes. 

I!.:}. :i  Two-Max wellian  Plasmas 

•  iarrett  showed  that  a  two-Max wellian  fit  is  often  a  good  representation  of  plasma 
distribution  functions  measured  during  geomagnetic  substorms. ^  The  density  and 
temperature  of  each  Maxwellian  component  can  be  found  from  four  velocity  moments  of 
the  measured  spectrum.  It  is  possible,  in  principle,  to  find  three-Maxwellian  fits  which 
match  sjx  moments,  although  the  effects  of  errors  in  measurement  of  the  plasma 
spectrum  become  increasingly  exaggerated  when  computing  the  high-order  moments.  It 
should  also  be  possible  to  find  multiplc-Maxwellian  least -square  fits  directly  from  the 
measured  distribution  'unctions  without  computing  the  velocity  moments  of  the  data. 


Singh'  He  mi  l.nergv 


V  t  .vo-Maxwelli  m  distribution  has  average  and  KMS  temperatures  given  by 


Vi  *■ 

ii  j  i  a. 
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(2-31) 


where  iij.  ii.,,  T|,  aad  T.,  are  the  respective  ileiisilies  and  temperatures  of  the  two 
components  of  the  spectrum. 

\  ingle  moiioi  nerget  ie  beam  can  match  two  velocity  moments  of  the  distributed 
spect rum  d  its  den  at  .aid  energy  are  chosen  according  to  l.qs.  (2- ltd  -  (2-22)  above,  for 
exainpl-  if  the  beam  density  is  equal  to  till'  total  plasma  density,  n.  *  n.,,  and  its  energ\ 
•s  I  y-'  then  the  zeroth  and  second  velocity  moments  of  the  two-Max  wellian  plasma 
and  the  n  niiociierget ie  beam  are  equal. 
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can  bo  chosen  individually  to  match  two  moments  of  each  of  the  Maxwellian  components 
of  the  spectrum.  In  this  ease*,  Fqs.  (2 -  lit)  -  (2-22)  would  tie  em|  loyed  along  with  the 
densities  and  temperatures  of  the  two-Max wellian  fit. 


The  second  approach  is  to  use  the  average  and  ItMS  temperatures  of  the 
two-Muxwollian  fit.  F.qs.  )  and  (2-31),  and  to  calculate  th<  beam  velocities  and 

densities  from  Fqs.  (2-2-1)  and  (2-2.r>).  In  both  cases,  as  many  as  four  moments  of  the 
two-Max  wellian  distribution  function  can  be  matched  by  two  mot  oonergetie  beams.  In 
practical  situations,  physical  considerations  would  be  required  to  make  a  choice  between 
the  two  methods  of  matching  velocity  moments. 


Three  or  More  Hearn  Energies 

The  moment-  of  a  two-Maxwellian  distribution  function  can  be  matched  in  several 
different  combinations  with  multiple  monoonergetie  beams.  As  it  the  two-beam  ease, 
each  Maxwellian  component  of  the  plasma  can  have  one  or  mor  beams  assigned  to  it 
which  individually  match  velocity  moments.  For  six-moment  :  latching,  three  beam, 
energies  and  densities  could  tie  selected  using  Fq.  (2-211)  for  each  component,  and  a  total 
of  six  beam  energies  would  be  required  to  simulate  the  two-Maxwellian  plasma.  As 
mentioned  above,  fir  computed  values  of  the  zeroth  through  fifth  moment  ol  the  full 
spectrum  can  aho  be  used  directly  to  find  three  beam  energies  and  densities  through  the 
iterative  minimization  procedure  described  in  Appendix  H. 


2.3.4  Arbitrarily  Assigned  Hearn  f.nergies 

The  velocity  moments  of  a  measured  distribution  function  can  also  be  matched  by 
monoonergetie  beams  whose  velocities  are  chosen  arbitrarily.  As  an  example,  four 
monoeaerget 'e  beam-'  can  match  four  velocity  moment  : 
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When  the  four  beam  velocities  are  fixed,  then  it  is  only  a  matter  of  solving  a  set 
of  linear  simultaneous  equations  for  the  beam  densities,  iij  through  n^.  It  should  be 
pointed  nit  that  not  all  combinations  of  bean)  velocity  may  be  chosen  for  v^  through  v^, 
because  negative,  and  therefore  unphysical,  solutions  for  the  beam  densities  cun  be 
obtainee  in  some  eases. 

T  ible  2-2  gives  the  densities  calculated  for  three,  four,  ana  i';ve  beams  as  a 
funetion  of  preassigned  beam  energies.  The  beam  energies  and  densities  are  normalized 
to  the  temperature  and  density  of  a  Maxwellian  distribution,  .and  the  velocity  moments 
used  for  the  calculations  are  given  by  the  right-hand  side  of  Kqs.  (2-8)  -  (2-1 2).  The  first 
three-beam  solution  in  Table  2-2  is  a  check  of  the  six-moment  solution,  Eq.  (2-29),  found 
by  the  iterative  prneed  ire  discussed  in  Section  2.3.2. 

I  ible  2-2  indie  des  that,  although  the  velocity  moments  of  the  monoenergetie 
ben  ms  are  matched,  tl  e  spectral  shape  of  the  beam  solutions  is  generally  not  similar  to 
that  of  i  Maxwellian.  For  a  simulation  facility  one  would  intuitively  prefer  an  envelope 
of  the  beam  density  which  roughly  approximates  the  Maxwellian  distribution,  Eq.  2-1. 
The  low  and  negative  values  of  density  found  for  some  combinations  of  beam  energy 
npparen'lv  result  from  forcing  the  beam  densities  alone  to  tiring  about  the  match 
between  the  velocity  Moments.  The  unphysical  and  intuitively  unsatisfying  results  using 
arbitral1  Iv  assigned  be  nil  energies  east  doubt  on  the  usefulness  of  this  approach  to  match 
velueit moments  of  d>  dributed  spectra. 


TABI.i:  2-2. 

PARTIChK  DKNSITIFS  FOR  PRFASSB 
BHAM  KNKRfilf-S 

;ni:d 

No.  of 

Beams 

Kj/kT 

n, 7 n 

0.303 

0.325 

3 

1.G57 

0.588 

4.931 

0.087 

0.5 

0.521 

3 

2.0 

0.339 

4.0 

0.140 

0.4 

0.282 

4 

0.8 

0.333 

1.f> 

0.009 

3.2 

0.316 

0.5 

0.518 

1.5 

0.150 

4 

2.5 

0.134 

3.5 

0.192 

0.2 

0.198 

0.0 

0.001 

5 

1.3 

0.548 

3.0 

0.127 

5.0 

0.000 

(1.2 

0.179 

O.fi 

0.102 

5 

1.5 

0.502 

3.0 

0.080 

5.0 

0.071 
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suction  :j 

S|>.\ri;cUAI'T  (MIMUMNC.  CALCULA  TIONS 


:i.i  K‘  \  i  ion  ai.i: 

111.-  previous  section  presented  some  mathematical  techniques  to  relate  the 
eharaeieri-dies  .if  undi  turbed  plasmas  to  those  of  one  or  more  monoencrgetic  beams  of 
charged  :>.irtie|es.  It  v  as  assumed  that  the  plasma  or  beams  proriuecd  a  flux  of  particles 
at  a  given  surface,  although  no  interactions  between  the  particles  and  the  surface  were 
console  red. 

In  this  -.notion  .ve  shall  compare  the  electrostatic  charging  produced  by  plasmas 
and  various  eombinatic  ns  of  monoencrgetic  electron  and  ion  beams  using  a  model  which 
accounts  for  several  of  the  interactions  between  the  incident  charged  particles  and  a 
"typical”  -paeecraft.  I' he  spacecraft  charging  calculations  arc,  of  course,  only  one  of 
several  possible  approaches  for  making  a  qualitative  comparison  of  the  effects  of 
plasmas  and  combinat  ons  of  monoencrgetic  brains.  A  spacecraft  simulation  facility, 
however,  will  devote  :  considerable  amount  of  its  effort  to  the  study  of  the  effects  of 
electrostatic  charging,  and  this  choice  for  comparison  can  be  justified  on  these  grounds. 

I!.:’  i'll  MtCIN'i;  MOD  X 

(i)) 

I'e  spaeeeralt  eharging  model  developed  by  (inrrett  calculates  the  equilibrium 
p"t<  1 0  ill  of  a  surface  which  receives  isotropic  fluxes  of  electrons  and  ions  with  arbitrary 

cm  u.  . . Ira  and  which  loses  charge  by  secondary  electron  emission,  electron 

back -cati'T,  and  photoelectric  ('mission.  The  model  has  been  rather  successful  in 
pfs’i-'lu  ;  the  potential  of  high- altitude  satellites  instrumented  to  measure  the 
di iTcfeut  >,  1 1  eni>r  ry  spe  -tra  of  electrons  and  protons  in  g<'o magnet ic  substorm  plasmas^ '  . 

Die  model  as  nines  that  the  spacecraft  can  be  represented  as  a  spherical 
l.nir'mair  ,u'o!><-  ,n  a  plasma  whose  I  leUve  length  is  much  greater  than  the  dimensions  of 
the  |»ro!.>.  Its  ei 1 1 t  ,  spectra  o|  the  plasma  electrons  and  ions  are  divided  into  l>'2 
<  m-rgv  '".ins",  out  tin  flux  of  charged  particles  to  the  surface  calculated,  taking  into 
account  the  ••tectro  tatic  potential  of  the  satellite  and  conservation  of  mass. 
M  i\ ,'.e| h  ((!,  two  Max  •.ellian,  and  arbitrary  spectra  observed  from  the  spacecraft's 
in  ■  t  rum/  a  I  it  i<  >n  can  !  x  h  »ulei !  into  the  energy  Ians. 


?0 


i 


i 


Secondary  electron  emission  from  electron  and  ion  bombardment  and  electron 
backscatter  are  calculated  as  a  function  of  the  incident  particle  flux  and  the  measured 
energy  dependence  of  the  secondary  emission  and  backscatter  coefficients  of  aluminum. 
Corrections  for  the  heterogeneous  surface  of  an  actual  spacecraft  are  made  by  small 
adjustments  of  these  coefficients  to  bring  the  calculated  potential  of  the  satellite  equal 
to  its  measured  value  when  the  satellite  is  in  "typical"  plasma  conditions.  Charge  losses 
by  photoemission  are  included  by  an  empirical  formula. 

We  have  modified  the  model  in  two  ways,  hirst,  the  time  dependence  of  charging 
was  included  by  representing  the  satellite  as  an  isolated  spherical  capacitor.  The  amount 
of  charge  gained  and  lost  by  the  surface  is  calculated  for  short  increments  of  time  in 
which  the  potential  is  held  constant.  The  not  gain  of  charge  is  then  used  to  compute  the 
new  value  of  potential  to  be  used  during  the  following  time  increment.  This  procedure  is 
repeated  until  the  potential  of  the  model  satellite  does  not  vary  in  succeeding  increments 
of  time. 

The  second  modification  was  used  only  for  potential  calculations  of  the  model 
when  irradiated  by  monoenergetie,  initially  parallel  beams  of  noninteraeting  charged 
particles.  It  accounts  for  the  electrostatic  deflection  of  the  beams  in  the  electric  field 
of  the  charged  body  which  attracts  oppositely  charged  particles  and  repels  particles  of 
the  same  sign. 

The  total  current  to  a  surface  of  arbitrary  shape  in  a  parallel  beam  is  simply  the 
product  of  the  current  density,  j,  and  the  geometric  cross  section,  A,  in  a  plane 
perpendicular  to  the  current  density  vector.  If  the  initially  parallel  beam  is  deflected  by 
a  symmetrical  potential  well,  the  deflection  can  be  represented  ns  an  "effective" 
cross-sectional  area  which  depends  on  the  strength  of  the  field  and  the  kinetic  energy 
and  charge  of  the  particles.  It  is  shown  in  Appendix  C  that  the  effective  areu  of  a 
spherical  conductor  of  radius  R  is, 

(m  <p  <  |qfM> 

(3-1) 

(qy?.  '  iqi-;|) 


t 


whore  \p  is  the  (si j^noti)  potential  of  tlie  sphere,  und  q  and  E  are  the  (signed)  charge  and 
initial  kinetic  energy  of  the  incident  charged  particles. 

For  the  charging  calculations,  the  electron  and  ion  current  to  the  model  satellite 
was  set  equal  to  the  sum  of  the  currents  from  the  monoenergetic  beams  each  of  which 
was  given  hv 


i  =  Ji  Aoft  (3-2) 

where  j.  is  the  unperturbed  current  density  of  the  i^  beam  with  energy  E.. 

The  secondary  emission  current  from  electron  and  ion  bombardment  and  the 
electron  backs  eat  taring  were  calculated  as  a  function  of  the  energy  of  the  incident 
particles  by  the  same  subroutines  used  by  (larrett's  model  for  distributed  energy  spectra. 

No  photoemission  was  included  in  the  spacecraft  charging  calculations  in  order  to 
simplify  comparison  of  the  results  between  monoenergetic  beams  and  distributed  spectra. 

3.3  It  KS  LETS 

The  spacecraft  charging  model  was  used  to  calculate  the  potential  of  a  spherical 
satellite  with  a  radius  of  1  meter  and  initial  potential  of  zero.  The  charging  by  plasmas 
with  several  different  electron  and  ion  temperatures  were  compared  to  charging  by 
beams  whose  energies  and  current  densities  were  selected  by  the  methods  discussed  in 
beet  inn  1.  Table  ;i-i  presents  the  parameters  of  some  of  the  Maxwellian  plasmas  and 
beams  and  for  the  charging  calculations. 


3.3.1  Single  Monoenergetic  beams  and  Maxwellian  Plasmas 

Charging  by  monoenergetic  electron  and  proton  beams  and  Maxwellian  plasma  was 
comput'd  for  several  beam  energies  and  plasma  temperatures.  The  current  densities  and 
energies  were  selected  so  that  the  first  (number  flux)  and  third  (enerjty  flux)  vcloeitv 
moment*  of  the  monoenergetic  beams  matched  those  of  the  Maxwellian  plasmas, 
Equations  (2  .’1)  and  (2-22).  For  this  ease,  the  beam  energies  were  twice  the 
corresponding  plasma  temperature. 
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(3-3) 
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TABU-:  3-1.  SPACECRAFT  ClIARCINC,  BY  MAX WELLIAN 
PI, ASM  AS  AN!)  MONO  BN  KR(I  KTIC  BBAMS 


PLASMA 


Lloctfons: 

Ions: 

Electrons: 

Ions: 


T  =  10  keV,  j  1.1)  iiA/rm2 
T.  =  10  keV,  j.  -  0.023  nA/eni 

I  =  10  keV,  j  =  1.0  nA/em2 
T.  =  20  keV,  j.  -  (1.033  nA/cm2 

l  Ji 


V?  -  -12.3  kV 

•q 


p  -  -14.2  kV 
<’(] 


BEAMS 


1  Electron:  I  =  20  keV 

■  G  a  .  2 

j  _  =  nA/em 

1  Ion:  L.  =  20  keV 

1  2 
j  =  0.023  nA/cm 

•  i 


1  Electron:  1'  =  20  keV 

e  2 
j  (  =  1.0  nA/cm 

1  Ion:  !..  =  40  keV 

*  o 

j.  -  0.033  nA/cm" 

Ji 


2  Electron’ 


2  Ion: 


el 

•*e2 
1  o2 
\>2 


=  5.69  keV 
=  0.41  nA/cm  2 
=  30.1  keV 
=  0.59  nA/cm2 


I...  =  5.69  keV 

1  ^  2 
jjl  =  0.0096  nA/cm 

I'...,  =  30.1  keV 
l2  9 

j  -  0.0 1 4  nA/cm 

Ji2 


p  =  -13.5  kV 
cq 


p  -  -1 2.9  kV 

I'M 


p  -  -15.0  kV 

cq 


t 
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\BL1-:  3-1.  (<  ’oncliidoii) 


3  Kloetron:  K  .  =  3.03  koV 

ot  , 

j  ,  -  0.16  nA/cin 
Jc  l 

K  „  =  10.6  kcV 
o2  2 

j^0  =  0.67  n  A  A*  m 

I'.  -  49.6  kcV 

oJ 

Jo3  "  n‘,T  V?  =  -11.9  kV 

3  Ion:  I',- ,  -  3.03  kcV 

"  2 
jj  -  0.0037  n  A  Amu 

.  10.57  kov 

j.  ,  0.0  16  n.A  'orn" 

Ji  1 

49.0  kcV 

I  •’  9 

j...  -  0.00  10  n A/cm' 

1 3 
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whore  n  and  T  are  t ho  density  and  temperature  of  the  Maxwellian  plasma  component,  q 
and  m  are  the  charge  and  mass  of  the  plasma  and  beam  parlicl  \s  (assumed  the  same 
species  4,  and  F^  and  arc  the  undisturbed  energy  and  current  <lensi  y  of  the  Imam. 

Figure  3-1  shows  the  charging  of  the  model  .satellite  wit  i  a  radius  of  I  meter 
under  irradiation  by  single  30  koV  electron  and  proton  beams  and  b  /  a  hydrogen  plasma  in 
which  the  electron  and  ion  temperatures  are  it)  keV.  It  can  be  ,c«*n  that  the  charging 
rate  and  equilibrium  potential  of  the  satellite  is  higher  when  exposed  to  the 
monocnergetie  beams,  although  some  differences  are  to  bo  ex[  acted  because  of  the 
important  influence  of  the  secondary  electron  emission  coeffie  cuts  on  the  ‘barging 
process. 

The  equilibrium  potentials  found  from  calculations  of  charging  by  Maxwellian 
plasma  and  beams  with  energies  and  current  densities  given  by  liquation  (3-3)  are 
compared  in  Figure  3-2.  The  correspondence  is  surprisingly  good,  considering  the 
crudeness  of  simulating  a  Maxwellian  velocity  distribution  by  a  single  monocnergetie 
beam. 

Figure  3-3  shows  calculations  of  charging  by  a  Maxwellian  plasma  with  an  electron 
temperature  of  10  keV  and  an  ion  (proton)  temperature  of  20  keV.  Charging  by  electron 
beams  with  an  energy  of  20  keV  and  proton  beams  of  40  koV  and  current  densities  for 
each  component  given  by  F.quation  (3-3)  are  also  shown.  In  this  ease,  the  eouilibrium 
potential  in  the  Maxwellian  plasma  is  somewhat  higher  than  under  irradiation  by  the 
beams. 

3.3.2  Two  Monocnergetie  beams  and  Maxwellian  l’lasma 

The  energies  and  densities  required  for  two  beams  to  match  four  velocity 
moments  of  a  Maxwellian  plasma  are  given  in  Fquation.s  (2-2?)  and  (2-28).  We  have 
calculated  the  charging  by  two  electron  and  proton  beams  and  in  Maxwellian  plasmas. 

Figure  3-4  shows  the  results  of  the'  calculations  for  electron  and  ion  beams  with 
energies  of  3.09  keV  and  30.1  keV  and  for  a  Maxwellian  plasma  with  electron  and  ion 
temperatures  of  10  keV.  The  equilibrium  potential  of  the  sa tel Ii t •  ■  model  is  more  than  2 
kV  greater  for  charging  by  the  beams  than  by  the  plasma,  although  the  charging  rate  is 
about  equal  for  both  cases  for  0  to  0.05  seconds. 
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SPACECRAFT  CHARGING 
BY  MAXWELLIAN  PLASMAS 
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FTGURK  3-4.  C'HARGTNC;  liV  MAXWF.I.I.IAN  l'T.ASMA 
AM)  TWO  KI.I.OTItON  AND  K'N  lil’.AMS 


Three  Monoenergelie  Mourns  and  Maxwellian  Plasma 

Charging  by  three  monoenerget io  eleotron  amt  three  moiiocncrge;  ic  electron  ami 
ion  bean  s  whose  velocity  moments  match  six  moments  ol  a  Ma.xwcl  aan  plasma  was 
computer,  usi  rise  the  spacecraft  charging  model.  The  beam  energies  amt  eurr  >nts  ware 
found  f r  ;n  liquations  (’’-29)  to  match  the  velocity  moments  of  a  M  iv.veihan  hydro, p-n 
pla-mia  vs  ah  an  electron  aim  ion  temperature  of  10  keV. 

The  results  of  the  charging  calculations  are  shown  in  Figure  There  is  very 

close  a. ."•cement  hetween  the  charging  rates  and  equilibrium  potentials  for  both  the 
throe-be  on  ami  Maxwellian  plasma  eases. 

1  F  *ur  li-niiis  am!  Maxwellian  Plasma 

T  e  char  fing  of  the  satellite  model  was  calculated  using  beams  chosen  to  •  imulate 
the  «iiff.  :it i.i !  ('nergv  spectrum  of  the  current  density  of  a  Maxwellian  plasma.  As 
discusser  in  Section  ’.2,  the  energy  distribution  was  broken  into  four  parts  and  the 
cl  rrent  •  ien-ulv  and  average  energy  of  each  part  computed,  using  l.quations  (2-4)  and 
(2-9). 

figure  shows  the  charging  using  the  four-beam  solution  given  in  Table  2-1 
compared  with  elnrgi”g  by  a  Maxwellian  plasma  with  electron  and  ion  temperatures  of 
]il  f  ,>\.  It  m  s« > m#' w  1  ta t  surprising  that  the  equilibrium  potential  l'oum:  with  four  electron 
mid  ion  'oimi--  chosen  to  mimic  tho  spectral  shape  of  the  Maxwellian  plasma  is  not  as 
close  as  with  other  en  es  with  fewer  licmns. 

H.'t.a  I •.•••> ?n ,  and  Two-Maxwellian  Plasma 

A  ,|j .cu.-  c,!  in  Section  9. .1.9.  the  velocity  distribution  of  a  non- Maxwellian  plasma 
cun  be  epproximated  I iv  a  t  wo -Max  welliau  distribution  function,  each  component  of  the 
distribution  being  ehurneleri/.cd  by  a  temperature  and  a  particle  density.  He  have 
e>  'input  *  ' !  the  charging  of  l  lie  satellite  model  in  a  plasma  with  a  t  wo-Max  wel  lian  electron 
(list ritnit loti  (unction  und  -inglo-Max welliau  ions.  I  ho  two  electron  eompouents  have 
t e 1 1 > : >■  r a 1 1 as",  of  1 1)  hoV  and  lit)  ke\ ,  ami  densities  ol  .1.0  ("in  and  0.4.1  cm  , 
respectively.  I’h* •  proton  plasma  has  a  temperature  of  10  keV  and  has  a  number  density 
equal  to  the  total  electron  density. 

We  have  compared  the  charging  by  the  l  wo-Max  welliau  plasma  to  that  of  several 
combinations  of  m* >m'oiiergot je  beams.  Table  .1-2  shows  the  beam  energies,  current 
densities,  and  resultant  equilibrium  potential  of  the  satellite  model. 


-.’SO- 


SPACECRAFT  POTENTIAL  (k’v) 


TIME  (seconds) 


FIGUKK  3-5.  CHARGING  HY  MAXVV  II,  1,1  AN  IM.ASMA  AND 
Til  K  KK  I : I . I •: <  ’TKGN  AND  ION  HLAMS 
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_  Maxwellian 

V?  = 

eq 

-12.5  kV 

_  Beams 

V  = 

eq 

-10.5  kV 

lOkeV 

TIME  (seconds) 


» if i i ! u i:  :t  r,.  '.iiauoinc,  nv  maxwkluan  im.asma  and 

I'fD'i:  I  | .  j;  ■!  I{t  )N  AND  ION  UK  A  MS 
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TABU-:  3-2.  SPArECRAlT  < 'I  I A  ll<  !IN< BY  TVVO-M  AXWELLIAN 
PLASMA  AND  MONOI'.N KIK I'/I  K 1  BEAMS 


PLASMA 


Electrons: 

1  -  10  keV,  j  .  -  0.8  tiA/cm" 

el  el 

T  0  =  30  keV,  j  =  0.2  nA/eni 

(OaV  =  «2.52«kiv 
(  Vrms  =  l4-°  kpv 

T.  =  10  keV,  j.  -  0.021  nA/cni 
l  Ji 

ri'<\ 

-1  1.0  kV 

Ions: 

BEAMS 

1  Electron 

B«  *  *Vh»w  •  “•»  k“v 

V 

eq 

=  -26.8  kV 

1  Electron 

l'e  -  3/2  (Te)  .w  =  18.8  keV 

1  Ion 

j  =  1.1  nA/cm" 

E.  =  3/2  T.  =  15  keV 

1  1  2 
j.  =  0.023  nA/cm 

eq 

=  -12.7  kV 

1  F.lectron 

l!e  =  2(Te»l!MS  !  28'°  k“V 

1  Ion 

j  -  1.0  nA/cm L 

Ej  =  2  Tj  =  20  keV 
jj  =  0.021  nA/cm2 

M'l  } 

-  -19.8  kV 

2  Electron 

E  ,  =  2  T  .  =  20  keV 
cl  el 

j  ^  0.8  nA/eni 

E  „  =  2  T  „  -  60  keV 

j  %  0.2  nA/cm2 
c2 

'P 

C<1 

--  -18.2  kV 

1  Ion 

R.  =  2  T.  --  20  koV 

1  1  9 

=  0.021  nA/cm" 

2  Electron 

Epj  =  7.92  keV,  =  0.54  nA/cm2 

E  „  =  51.9  keV,  j  0.46  nA/cm2 

c2  •’cz 

=  -19.2  kV 

2  Ion 

E.,  =  5.69  keV,  j..  =  0.0096  nA/cm 

il  Jil  7 

Ej2  =  30.07  keV,  jj2  -  0.014  nA/cm" 

Cq 
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li . quibbi-mm  potential  found  with  a  single  electron  beam  i>.  presented  >n  show 

the  effect  nl'  •  emovmg  ions  from  tin-  simulation.  Without  the  i<  n  component,  the 
satellite  mi. del  charges  until  the  secondary  electron  emission  and  back  scat  ter  are  equal 
to  the  iiiMd-  nt  electron  flux.  The  equilibrium  potential  is  close  to  t  int  of  the  electron 
beam  because  the  secondary  electron  emission  coefficient  peaks  at  m  iticrgy  of  a  few 

(M) 

hundred  i  h-etronvolt  ■  '  and  is  small  at  higher  energies. 

The  sm^le-oleetron  and  single- ion  beam  energies  and  currents  m  ,  uo.e  3-2  .vi-r. 
chosen  to  match  two  velocity  moments  of  the  iwu-Maxweiiian  plasma,  as  discussed  in 
Section  '.bit.?.  The  two-electron  and  single-ion  beam  energies  and  currents  match  the 
fir.t  and  third  velocity  moments  (particle  and  energy  flux)  of  each  component  of  t he 
distrihutic.il  functions. 

The  energies  and  currents  of  the  tivo-cjretroii  and  two-ion  beam  case  were  found, 
using  lariat  ions  (2-2-1)  and  (2-?a),  to  match  four  velocity  moments  of  the  distribution 
function-.  ba.-ed  on  the  average  and  It. MS  temperatures  of  the  plasma  particles. 

Til--  discrepancies  between  the  calculations  of  equilibrium  potential  in  the 
two-Maw.ellinn  plasma  and  in  moiioenerget ie  beams  are  somewhat  greater  than  tho-e 
found  with  a  single  Maxwellian  plasma.  I'be  difference  limy  be  caused  by  the  higher 
temperature  component  of  the  electron  plasma,  which  skews  the  second  and  third 
velocity  moments  of  the  electron  distribution  function.  The  high-energy  electron  beams 
require.!  to  match  these  velocity  moments  apparently  have  a  strong  influence  on  the 
equilibrium  potential  of  the  model. 

.’!.t  i  >1  S<  I’SSION 

The  calculation-,  give  a  qii.aiit.ative  idea  of  the  charging  which  would  be  observed 
in  a  ”.•>:•  m •era ft  testfg  faeilitv  in  which  mortoenerget ic  beams  were  used  to  simulate 
space  plasmas  with  distributed  energy  spectra.  As  expected,  the  equilibrium  potential  « >1 
the  spacecraft  under  test,  and  therefore  the  charge  density  on  its  surface,  is  only  a 
function  if  the  electron  and  ion  beam  energies  and  currents.  An  important  result, 
however,  U  tin'  observation  that  the  inonoenergetie  beams  can  be  chosen  to  match 
several  velocity  mono  nts  of  a  distributed  spectrum  and,  at  the  sunn  time,  produce  the 
same  china'  density  on  the  spacecraft.  Tims,  surface  phenomena  which  are  influenced, 
for  example,  by  energy  flux  as  well  us  charge  density  can  be  investigated  m  a  laboratory 
facility  with  -i  reasonable  degree  of  confidence  m  the  simulation  fidelity. 


It  should  bo  made  clear  that  the  charging  model  list'd  here  is  a  very  simple  one  and 
does  not  account  for  the  complex  geometry  or  surface  details  of  a  real  spacecraft.  More 
complicated  ehargin  '  codes  exist,  however,  which  could  he  used  to  make  more  detailed 
comparisons  of  spacecraft  charging  by  monoenorgetie  beams  and  space  plasmas.  The 
NASCAP  eode^111,  t\^r  example,  is  probably  the  most  ambitious  attempt  to  represent  the 
geometrical  and  surface  configuration  of  real  satellites  in  the  environment  of 
geosynchronous  orbit.  Modifications  of  NASCAP  would  be  required  to  calculate  the 
charging  of  a  three-dimensional  object  tinder  irradiation  by  beams  of  charged 
particles/12'  but  it  is  likely  that  NASCAP  would  be  a  useful  tool  for  comparing  the 
conditions  of  laboratory  simulation  to  those  of  space. 


:t'>- 


SUCTION  1 

SUMMARY  AND  K ROOM M UN D ATlONS 


4.1  SUMMARY 

Wo  have  examined  matliemat  ical  techniques  to  choose*  the  .-nergy  and  currant 
density  of  monoencrgelie  beams  to  simulate'  the  distributed  spectra  of  p.asmas  in  space. 
In  the  first  approach),  the'  differential  current  density  spectrum  of  the  plasma  was  divided 
into  a  number  of  energy  bands.  The  beam  energy  and  current  were  calculated  for  each 
band  to  provide  a  piecewise  reproduction  of  the  distributed  spectrum.  This  technique  ls 
probably  tin'  most  intuitively  satisfying  when  a  large*  number  of  beams  with  different 
energies  can  be  used  for  the  simulation  because  the  envelope  of  the  velocity  distribution 
of  the  bourns  e an  <»lose)y  mimic*  that  of  the  plasma. 

The  second  general  approach  was  to  choose  the  beam  energies  and  current 
c b 'ii si t i •  -s  to  match  the  velocity  moments  of  the  plasma  distribution  function.  The* 
veloeity  moments  ore-  averages  related  to  physical  quantities  such  as  particle  density, 
flux,  pre  ssure,  and  energy  flux,  and  have’  been  used  extensively  to  characterize  the 
measured  properties  of  plasmas  in  space.  We  have  found  expressions  for  beam  energies 
and  densities  in  terms  of  the  plasma  properties  such  as  the  average  and  RMS 
'hi'inperd uro",  density,  and  average'  velocity.  Combinations  of  one',  two,  and  three’ 
beam  wer»*  found  to  mate’ll  two  to  six  veloeity  moments  of  Maxwellian  distributions. 
I  be  see ■■  t • ' « • ! u * i < nil's  also  can  tie*  applied  to  other  spectral  shapes,  and  they  wore  usee!  to 
ex;. mite  t  wo  Maxwellian  distributions.  Unphysiea!  or  intuitively  unsatisfying  result'. 
'■vMv  found  wb'-n  the’  problem  was  overspeeifie'd  by  arbitrarily  selecting  beam  energies 
and  b.  ;e u  deii-ite"  calculate!  I  to  mate’ll  a  set  of  velocity  moments  of  a  distributed 
s;*ee ‘  t  I’ll : a . 

A  rapf-  computational  model  was  iise'e!  to  e’onipare  the'  charging  of  a  space'Cl’a  ft 
tiv  plauicr  with  distributed  spectra  and  by  uioneie'iiergetio  beams.  The  plasmas  were 
siuular  to  tlio  <’  toimd  m  high  orbits  (luring  gi'omngnct ic  substoniis.  with  multikilovolt 
b  mperat nr"s  and  a  I'.-w  partud"  .  per  cubic  cent  iinet er.  These'  calculations  were  made  to 
gain  a  qualitative  comparison  of  the  approaches  for  choosing  monoe'iierge't ie  beams  to 
'  miniat"  sp  ■.•(*  plasmas. 

\  1 1 1 1  .iigt.  a  e|o.o  eom.i.iri  on  was  not  . . .  when  only  a  few  beams  were  used 

t"  •iiimlnt"  the  di  I fat.u I  "<i  ./ i.  ei rum  ot  a  pla  mu,  sonic  combinations  of  beams  gave 
’imilar  "his'i'ing  rat"  and  e  qiiibi.riiim  potential  ..  The  equilibrium  potentials  found  using 
b"aai  to  male’ll  v.looiK  moan  ats  of  a  t  wo- Ma  x  w  el  halt  plasma  generally 


were  within  ;i  few  kilovolts  ol  charging  l*y  the  distributed  spectrum,  but  showed  more 
divergence  than  the  simulations  of  simple  Maxwellian  plasmas.  The  spacecraft  charging 
model  cannot  be  teed  as  the  only  criterion  to  evaluate  the  fidelity  of  simulation  of  a 
given  plasma  environment,  because'  physical  considerations  other  than  the  charging  rate 
and  equilibrium  pot'  ntial  must  also  be  considered  in  this  evaluation. 

4.2  Rl.r OM.MKNT)  \T!ONS 

A  number  of  characteristics  of  space  plasma;;  have  not  been  included  in  this 
analysis  which  may  be  important  for  certain  simulation  situations.  The  most  evident  of 
these  is  the  observed  anisotropy  of  the  velocity  distributions  of  particle's  in  space. 
Plasma  drifts  are  common  as  well  as  distributions  which  show  strong  correlations  with 
the  direction  of  the  magnetic  field  in  space.  Plasmas  are  often  characterized  as 
"drifting  Maxwellian",  "bi-Maxwellian"  (temperature  parallel  to  the  magnetic  field 
different  from  perpendicular),  or  "loss  cone"  distributions  (very  few  particles  in  a  cone  in 
velocity  space).  These  anisotropies  can  be  important  for  interactions  with  spacecraft 
with  large  oriented  surfaces  or  cavities  aligned  with  the  directed  plasma  flux. 

We  recommend  that  the  present  work  be  extended  to  account  for  anisotropies  in 
the  velocity  distributions  in  space.  A  number  of  techniques  exist  to  describe 
mathematically  the  anisotropies,  such  as  expressing  the  distribution  function  in  terms  of 
spherical  harmonies.  Some  of  these  techniques  can  probably  be  used  to  provide  objective 
relationships  between  plasma  parameters  and  the  placement  of  electron  and  ion  sources 
in  a  simulation  facility. 

Except  for  the  calculations  of  average  charging  using  the  simple  spacecraft 
charging  model,  we  have  not  considered  the  interactions  between  the  plasma  particles 
and  the  surface  of  a  spacecraft.  In  fact,  the  behavior  of  a  dielectric  material  irradiated 
by  electrons  and  ions  depends  to  a  large  degree  on  the  profile  of  charged  particles 
trapped  in  the  upper  atomic  layers  of  the  surface.  This  profile  is  a  function  of  the 
material,  particle  species,  energy  distribution,  mid  angle  of  incidence  on  the  surface.  An 
isotropic  plasma  with  a  distributed  energy  spectrum  will  produce  a  very  different 
charging  profile  from  that  of  a  monoencrgctic  beam. 


Ur  roeomiiiend  lint  further  theoretical  and  computational  studios  lie  performed 
to  il  lir.l  late  these  di  I  I  <  rouees  and  provide  guide  lilies  (Or  the  selection  ol  uionoener  got  ie 
Sienins  to  reproduce  the  charging  effects  from  plasmas  in  typical  spacecraft  material-  .  \ 

mimher  af  computer  codes  exist  which,  through  Monte  Carlo  routines,  follow  the 
trajectories  of  "typical"  energetic  particles  in  materials.  They  account  for  deflections 
by  elast  e  and  inelastic  collisions  and  energy  losses  by  the  same  processes  and  b\ 
radiation.  1'he  equilibrium  profile  of  charge  in  art  insulating  material  can  be  found  oy 
following  the  trajectories  of  a  few  thousand  particles  and  simple  models  of  chargee 
particle  nobility  in  the  material.  The  eodes  can  accept  a  variety  of  input  conditions, 
including  distributions  of  particle  velocities  and  incident  angles,  and  can  be  run  on  a 
moderate  size  computer. 

Calculations  of  the  interactions  of  charged  particles  with  spacecraft  materials 
should  be  coupled  with  experimentation  to  test  their  validity  and  verify  results  and 
prediction-.  We  recommend  that  a  small-scale  experimental  program  be  undertaken  to 
accomplish  these  ends.  Samples  of  spacecraft  material  should  be  simultaneous  irradiated 
by  electrons  of  more  than  one  energy,  ions  (preferably  protons)  and  light  with  a  strong 
eompotu  at  a  the  extreme  ultraviolet.  The  experiments  should  be  performed  under  very 
good  vacuum  conditions  to  avoid  surface  contamination  by  oil  or  ol  or  materials  which, 
are  atypical  of  the  space  environment. 


Seim-  exp'Timentation  of  this  type  has  already  been  reported,  by  the  spacecraft 

(e  o  :}) 

charging  -nuontilie  community.  ”  llowevr,  the  materials  were  usually  irradiated 
by  a  single  electron  beam,  and  no  attempt  was  made  to  correlate  surface  phenomena 
with  theoretical  analysis.  The  fact  that  no  ion  beams  were  included  in  the  tests  may  also 
strongly  modify  the  surface  interactions  in  comparison  with  those  in  the  plasma 


environment  of  space. 
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APPENDIX  1! 

three-beam  solution  to  match  six  velocity  m<  >m i ■; n ■  is 


To  find  a  three-ben rn  solution  to  match  six  velocity  mom- fits  of  a  distributed 
spectrum,  six  nonlinear  simultaneous  equations  must  be  solved.  ')  !,  -si-  equation.,  are  of 
the  form 


£  £ 

k=l  j  0 


where  n,  and  v.  are  the  densities  and  velocities  of  the  beams,  and  M.  ropn»M-  its  the  j 
K  K  J 

velocity  moment  of  the  distributed  spectrum. 

Although  any  physically  reasonable  distribution  function  can  be  utilize, we  used 
the  six  velocity  moments  of  a  Maxwellian  distribution,  Lqs.  (2)-(7\  to  find  s  dutions  of 
Eq.  (B-l).  The  beam  densities  and  velocities  were  normalized  by  cliai  giruj  variat  les. 


ni,:v? 


v  i  ,  ,  i 


that  the  moment  equations  for  a  Maxwellian  distribution  become 


a  -I-  b  +  c  -  1 


.ax  l  by  I  cx 


2  .  .  2  ,  Z  _  In 
ax  I  by  •(  ez  - 


>  ']  J  7T 

ax  f  by  f-  ex  - 


,  A  ,  h  lr>  2 
ax  +  dv  +  c.v.  =  -  v  n 
(>'! 


rj  J  ,  fj  .  'j  3  2 

.lx  S  l  > v  4  cy.  -  ■  it 
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K<|:.  <!!-.'!)  can  |n>  solved  by  ;m  iterative  procedure  that  minimizes  the  expression 
represent",!  bv  the  absolute  sum  of  Kqs.  B-3.  Trial  solutions  are  substituted  into  u 
computer  program  which  then  converges  on  the  best  solution  through  an  iterative 

(r>- 1 ! 

process. 

The  only  solutions  of  Kps.  (li-3)  found  by  the  iterative  routine  were  permutations 
of  the  following: 


( l.K,OHr),  0.AK81) 
(a, !.,e)  (n.osoi.,  (i.  o.  3:::irj) 


CL5-4) 


The  energies  of  ‘he  beams  were  found  from 


(U-5) 


and  si  in  i  1  >r  expressions  for  1'.,.  and  H 

x-  «) 


I'.,!.  Bevmgton.  I >;jV i  liednelion  and  hrror  Analysis  for  the  Physical  Sciences 
Cbehraw -Hill  It,. ok  Company,  New  York,  l'JGil),  Ch.  11. 


APPENDIX  (' 

Fl'FHf'TIVi:  AREA  OF  Si'll  ERK  *AI,  niNDl'CTOH 


We  shall  assume  that  a  sphere  of  radius  R  is  iti  a  uniform,  parallel  bejtm  of 
non-interacting  charged  particles  of  mass,  m,  and  (signed)  chart  e,  <].  If  trie  sphere  is 
uncharged,  the  eurrent  to  the  sphere  is  given  by, 


1(0)  =  7T  R”j 
■  o 


(r-i ) 


where  i  is  the  eurrent  densitv  of  the  undisturbed  heron. 

o 

When  the  object  charges,  the  particle  trajectories  of  the  beam  are  deflected  by 
the  static  electric  field  around  the  sphere.  The  effective  area  of  the  sphere  is  then  the 
circular  area  whose  radius  is  the  "impact  parameter"  of  the  beam  particles  that  just 
graze  the  surface  of  the  sphere. 

Elementary  texts  of  elasMonl  mechanics  show  that  the  impact  parameter,  b,  is 
related  to  the  distance  of  closest  approach,  a,  according  to  the  relation, 

b2  + 

1  1  (C-2) 

a  —  irv 

<  j 

9 

The  initial  kinetic  energy  of  the  particles  is  l/2mv”  and  (a)  is  the  electrostatic 
potential  at  the  distance  of  closest  approach  to  the  sphere.  Setting  the  kiiutic  energy 
equal  to  |qV  |,  a  =  the  effecti  e  area  of  the  sphere  is  then, 


v  b 


7T  K'  ( t  -  — ~~  )  ('*  ;  <!',V„l) 

MVJ 


(('-.la) 


When  the  potential  of  the  sphere  is  greater  than  the  initial  aeeolorat ing  voltage  of 
like-signed  particles,  the  particles  do  not  reach  the  surface. 

A,.f  r  (m  y>.;  ■:  |-iv  J  )  (r-:ih) 
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